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ABSTRACT 



Context. Galaxy clusters are among the best targets for indirect dark matter detection in 7-rays, despite the large 
astrophysical background expected from these objects. Detection is now within reach of current observatories (Fermi- 
LAT or Cerenkov telescopes); however, assessing the origin of this signal might be difficult. 

Aims. We investigate whether the behaviour of the number of objects per 'flux' bin (log N — log F) and that of the 
stacked signal could be used as a signature of the dominant process at stake. 

Methods. We use the Clumpy code to integrate the signal from decaying or annihilating dark matter and from cosmic 
rays along the line of sight. We assume the standard Navarro-Frenk- White (NFW) profile for the dark matter density 
and rely on a parametrised emissivity for the cosmic-ray component. In this context, the consequences of stacking are 
explored using the MCXC meta-catalogue of galaxy clusters. 

Results. We find the value of the slope of the log N — log F power law (or the increase of the signal with the number of 
stacked objects) to be a clear diagnosis to disentangle decaying dark matter from cosmic-ray induced 7-rays. For dark 
matter annihilation, depending on the signal boost from the substructures, it is either similar to the cosmic-ray (CR) 
signal (no boost) or similar to the decay case (large boosts). The shift between the brightest object and its followers 
also depends on the signal origin. For annihilation, this shift and the stacked signal are poorly constrained because of 
the large uncertainty affecting the boost. We also underline that the angular dependence of the annihilation signal is 
not universal because of the substructure contribution. 
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1. Introduction 

Contamination by astrophysical backgrounds is the curse of 
dark matter (DM) indirect dete ction searches. Ways around 
this issue are mult i- wavelength (iRegis fc Ullioll2008l) or/and 
multi-messenger (|Pato et al.l |2009[ ) analyses. In 7-rays, 



in several of these objects: P erseus and Abell 2029 from 
Whipple (iPerkins et al.ll2006l) : Abeh 3667 and 4038 from 
Kangaroo lKiuchi et al.l I2009D: Coma Abell 496 and 85 
from H.E.S.S. (lAharonian et al.|[2009allb[ ). A broader sur- 



background-free objects su ch as dwarf sphero idal galaxies 
(dSphs) are another option (j Walker et al .1120111 ). For targets 
presenting an astrophysical background, a useful diagnosis 
comes from the spectral behaviour, which differs for DM- 
induced or CR-induced 7-rays: DM signal s should exhibi t 
a sharp cut-off at the DM mass (.Bringm ann et al.ll2011|) . 
whereas astrophysical processes produce power laws with a 
slope ~ 2 — 3 (the cut-off energy depends on the source). For 
instance, the Galactic centre signal obse rved by the High 
Ener gy Stereoscopic System H.E.S.S. (jAharonian et al.) 
120041) is consistent with astrophysical sources. Despite their 
astrophysical background, galaxy clusters are also promis- 
ing DM targets in 7-rays (or at radio and microwave fre - 
quencies) as first proposed in IColafrancesco et all (|2006[) . 
Cerenkov instruments have set constraints on a DM sig- 
nal or astrophysical background from the non-detection 
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vey has been made by the Fermi-LAT i nstrument, where 
33 (re sp. 49) clusters were analys ed in lAckermann et all 
(|201G[ ) (resp. lAndo fc Nagail [2013) . T he consequences for 
hadro nic models are discussed, e.g. in iJeltema fc Profumol 
()201lD . Recently, a detection was reported by the MAGIC 
collaboration in the Perseus Cluster. However, this signal, 
also detected by Fermi-LAT (jAbdo et al.ll2009l ). is consis- 
tent with an emission f rom the giant elliptical ga laxy NGC 
1275 lying at its cen tre (lAleksic et al.ll201(]ll2012[) . From the 
Fermi-LAT results. fColafrancesco et all ( 20101 ) argue that 
it should be possible to resolve and detect the diffuse 7-ray 
flux of astrophysical origin coming from t he outer co r ona o f 
the Perseus cluster. Even more recently, iHan et aP (J2012l) 
reported a detection for Virgo, Fornax, and Coma from an 
analysis of Fermi-LAT data. These authors find that an in- 
terpretation in terms of DM annihilation is preferred by 
their analysis, although a CR origin is also possible. 

From these promising results, we anticipate that many 
clusters will be detected with the current or the next gen- 
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eration of 7-ray instruments. Multi-wavelength analyses, 
search for a spectral feature, and angular dependence of 
the signal arc tests to disentangle 7-rays of astrophysical or 
DM origin. In this paper, we investigate another test based 
on the study of a population of sources, thus somehow po- 
sitioned between single source an d full sky power-spectrum 
analyses (jAndo fc Komatsull2006[ ) . Actually, several studies 
hav e compared the astrophysical and DM a nnihilation sig- 
nal (jJeltema et al.ll2009t (Pinzke et alll201lD as well as DM 



annihilation and decay (JHuang et al.ll2012l ) on a somewhat 
limited number of clusters (due to the size of the available 
catalogues at the time). However, the systematic study of 
the prospects of a stacking analysis with respect to single- 
source analyses has not been investigated before. 

This paper is part of an ongoing effort to address this 
question in the context of the recently as sembled MCXC 
meta-catalogue of 1743 X-ray clusters (jPiffaretti et al.l 
l201lh . A first paper in this series highlighted the poten- 
tial improvement (a factor > 5 — 100) that can be brought 
by a stacking a nalysis over a single- source analysis for the 
DM decay case (jCombet et al.ll2012[ ). A second paper con- 
sidered the DM annihilation case in the light of the sensi- 
tivity of existing and future 7-ray instruments (jNezri et al.l 
I2OI2I ). It showed that stacking only brings a factor of '^ 2 
improvement in sensitivity for the Fermi-LAT instrument, 
but no improvement for Chcrcnkov Telescope Array (CTA)- 
like observatories. This third and last paper of the series 
inspects the case of the CR signal to address the possi- 
bility of using a stacking strategy to disentangle it from 
DM-induced signals. For that purpose, we use a generic de- 
scription of the CR signal and refer the reader to the two 
previous papers (jCombet et al.ll20ia : iNezri et al]|2012D for 
a detailed description of the DM modelling. We stress that 
this study remains at the phenomenological level. We are 
waiting either for the detection of several objects or a bet- 
ter knowledge of the input ingredients in order to quantify 
further the potential of the new diagnosis we investigate. 

We present the ingredients and the calculation in 
Sect. [2 We then present the results in Sect. |3] and conclude 
in Sect. H 

2. Ingredients and calculation 

The total flux expected in a given direction (l, b) (Galactic 
coordinates) and integrated over the solid angle AD, is given 
by the product of an energy-dependent term with an astro- 
physical term A^ 



d(t){E,i,b,An) 

dE 



dN 

— {E)xA{l,b,An). (1) 

In this study, we discard the spectral term and focus on 
A, which encodes all the information about the spatial de- 
pendence of the signal, and the relative intensity between 
clusters. A will be termed 'flux' in the following, but should 
be understood as the astrophysical contribution to the ac- 
tual 7— ray flux. For an observation in the direction of the 
cluster's centre {Id, &ci) and using an integration angle a-mt, 

Ax{aint)=Ax{lchbch^^)^ f [ £x{l',^)dl'dn (2) 



^ The redshift distribution of the MCXC catalogue of galaxy 
clusters (Piffarctti ct al. 201 j) p eaks at z ^ 0.1 (s e e thei r Fig. 1). 
Following .Coinbct ct al.. (,20ia ) and iNezri et"aLl (|2012l ') , we ne- 
glect the absorption for the MCXC galaxy clusters, so that the 
energy-dependent term and the spatial term are decoupled. 



is the integral of the 'emissivity' £x{l' , ^), over line of sight 
I' and solid angle Ail = 27r-(l — cos(aint)) for the process X. 
For DM, ^Docay = RDM for decay and fAnnihii. = Pdm for 
annihilation, where poM is the DM density profllc in the 
galaxy cluster, as presented in Sect. 12.11 For CR, fcR = 
CcR, as defined in Sect. 12.21 

The integral Eq. ([2]) is comput ed for all MCXC clusters 
with the CLUMpMEfcode v2011.09 (|Charbonnier et al.l20ll) 
adapted to include the CR-induced 7-ray case described 
below. 

2.1. Dark matter halos 

The DM distribution for each galaxy cluster is parametrised 
from fsee lCombet et aDl2012l for more details): 

— its Msnn value ^ pro vided in the MCXC catalogue 



(jPiffaretti et al 



mM; 



the choice of a un iversal DM profile, here the w idely 
used NFW profile (jNavarro. Frenk fc WhitHl997t) . 



PNFwir) 



(3) 



where rs is the scale radius and ps is the normalisation; 
a mass-concentration relationship, where the concen- 
tration is defined to be ca = Ra/ts for an NFW 
profilqj. The la tter is observationa lly constrai ned at 
the c l uster scale (fPointecouteau et ah 2005; Buot e et al.l 
l2007t lEttori et al.ll2010t) . and it has al so been exten- 
sively studie d in numerical simulations (IWechsler et aLl 
2002, 2 0061: iZhao eTaD [2001 [200I iNeto et al.l 120071: 
Gao et airB oOSl. 120111: 'Maccio ct al."2008': 'Duff v et al ' 



2008, 



2011 



2OIOI: iGiocoh et ah 2010: Mvuioz-Cuartas et al 



Klvpin et al. l 20111). In this s tudy, we use the 



parametrisation of iDuffv et al.l ()2008[ ). 

For DM annihilation, we also have to take into ac- 
count the contribution from substructures as they have 
been shown to boost the signal. In this work, we assume the 
following for the mass and spatial distribution of the sub- 
structures: i) dNsuhs/dM (X Af ~"^' w ith qm = 1.9, a mas s 
fraction /dm = 10% in substructures (jSpringel et al.ll2008l ). 
a minimal and maxim al mass of 10~^ Mp) a nd 10~^Mciuster 
respectively, and the iBullock et al.l (|200ll ) concentration 
(down to the minimal mass); ii) the substructure spatial 
distribution dNsuhs/dV follows the host halo smooth pro- 
file. Recent high-resolution numerical simulations that are 
for for Galaxy-like objects (Springel et al. 2008i) or are dedi- 
cated to cluster-size halos (jGao et al.ll2012b[) provide differ- 
ent figures for the slope aM , the mass fraction /dm , and the 
spatial distribution of substructures. We will briefly com- 
ment on how other choices impact on our analysiqj. 



^ |http : //Ipse ■ in2p3 ■ f r/clumpy/j 

^ Ma (with A = 200 or 500) is the mass enclosed in a sphere of 
radius -Ra, the radius within which the average density reaches 
A times the critical density of the Universe. 

^ In practice, w e work with virial quantities (see, e.g. 
iGiocoIi et al.|[2010l . to switch from A to virial quantities), e.g. 
Rvii is the physical size of the galaxy cluster halo. 

^ A detailed discussion of the annihilation signal and its 
dependence on the subs t ructu re parameters can be found in 
Sect. 3.4 of INezri et al.l (|2012l 'l and is not repeated here. We 
underline that taking the same distribution as the smooth halo 
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2.2. Cosmic-ray component 

High-energy 7-rays can also be produced from astrophysical 
processes, in particular from the inter action of cosmic-ray 
proton s with the gas of the cluster. iPinzke fc Pfrommerl 
(|2010f ) found a universal radial dependence of the emis- 
sivity fcR = C'cr('') of the clusters based on cosmologi- 
cal simulations. These authors acknowledge, however, that 
the CR spatial distribution from which 7-rays are emit- 
ted could be affected by, e.g. additional CRs injected from 
AGN, or CR diffusion in m omentum and space. For in- 
stance, lEnfilin et al.l (|2011[ ) argue that merging clusters 
should have a more centrally concentrated CR popula- 
tion than relaxed ones, providing bi-modality of their 7-ray 
emissivities. Nevertheless, we limit ourselves to the simplest 
case and assume a universal dependence across the cluster 
population we study here. We follow the fo r malism from 
IPinzke fc Pfrommeiil (|2010[ ) and lPinzke et all ()2011f ): 



CcR(r) 



C200 = 1.7-10 



' C200 - 5 ■ 10-7 
1 + (h±7) 



-6 



5-10" 



PgasW, (4) 



.7/ M, 



200 



ViOI^Mg 



0.51 



1.04 



M, 



200 



and fltrans = 0.021 R 



200 



M: 



lOi^Mp 



0.39 



0.15 



200 



1015M, 







Similar to the X-ray luminosity, the 7-ray emissiv- 
ity is proportional to the square of the gas density 
Pgas(^)- It is thereby very sensitive to the densest cen- 
tral parts of clusters. The relation between the DM and 
Pgasir) has been extensively studied from a theoreti- 
cal point of view (iMakino et al.l Il998t ISuto et al.l [TqqI 



It lOstriker 
[[[2009; Ca 



et al.l 120051: iBode et alj 



Komatsu & Scli akI 120011: 

2009; FrcderiksQiLeLalJ [2009; C apelo et al.l [2012D From 
the REXCESS representative sample. ICroston et al.l (|2008( ) 
showed the universal behaviour of Pgas{r) in clusters 
(though a large sca tter exists at the centre) . The AB model, 
first introduced bv lPratt fc ArnaudI (|2002f ). and given by 



/9gas(a;) OC 



-3fc2/2+fci/2 



(5) 



with X = r/i?5 oo, was fitted to the average REXCESS den- 
sity profile by lArnaud et al.l (|2010f ) , with best- fit parame- 
ters Xc = 0.303, fci = 0.525, and fc2 = 0.768. This analytical 
phenomenological model reproduces well the cuspy shape 
of density profiles derived from X-ra y observations and wa s 
also used for the MCXC assembly (|Piffaretti et al.ll20lH) . 
In order to remain consistent with respect to the data we 
use, the same AB profile is used here. The normalisation of 
Eq. dS]) is chosen to give fgasM^m at -^500- From the MCXC 
masses M500 and the empirical relatio n between /gas an d 
M500 given for the REXCESS sample ([Pratt et al.n2009D . 
an average value of /gas = 0.10 ± 0.06 is derived, which 
we use for the whole MCXC sample. We note that the gas 
fraction increases with radius in halos (to ultimately tends 
towards the cosmic value), and the intrinsic scatter is im- 
portant and scales with the cluster mass JVikhlin in et al.l 
I2OO6I: lEttori et al.ll2009l:IPratt et al.l [2001 1 201 Oh . 



for this quantity or using that of iGao et al.l (l2012bf) does not 
impact the annihilation results (see lNezri et aI.ll2012D . 



3. Results 

We first comment on the universality of the 7-ray signal 
angular dependence. In particular, we inspect whether this 
universality holds for any signal hypothesis (annihilating 
DM, decaying DM, or CRs); in that case, observations of 
different objects may be optimally stacked to increase the 
signal. We also inspect whether a given angular dependence 
can be associated to a unique signal origin. We then turn 
to the comparison of the galaxy cluster population signals. 



3.1. A different angular dependence for different origins? 

The spatial dependence of the signal has been suggested as 
a test to disenta ngle DM decay from annihilation, in par - 
ticular in dSphs (|Palomares-Ruiz fc Sieeal-Gaskinsl [20101 ). 
Indeed, the use of universal DM profiles ensures that the 
projected profiles are universal if properly normalised by 
their value at r ~ r.,. By construction, the projected pro- 
files for decay are universal, and so they are for annihilation 
if the smooth signal onl y is taken i nto acc ount. However, 
as can be de duced from [Gao et al.l ([2012^ results and as 
underlined in lNezri et al.l (|2012[ ). the substructure contribu- 
tion overtakes the smooth contribution at various distances 
from the centre, depending on the mass of the object rather 
than on its scale radius. Small objects (like dSphs) are dom- 
inated by their smooth contribution up to large distances 
from their centre, whereas massive galaxy clusters are al- 
most completely dominated by the substructure contribu- 
tion. As a result, departure from universality of the total 
DM annihilation signal (the one which we have access to) 
is expected. 

This is illustrated using all MCXC galaxy clusters. 
Defining for any cluster its virial angle avir = tan~"'^(i?vir/'^) 
and its scale angle Us = tan~^ (rs/d) (d its distance to the 
observer), we show in Fig. [T]the fraction of the integrated 
signal with respect to the maximal signal obtained, i.e. 



Ix{il) 



-4x(avir) 



(6) 



as a function of ry = ajnt/as- 

For DM annihilation, the smooth contribution (left- 
hand envelope) is the most centrally concentrated as it 
comes from the integration of the steep DM profile squared 

VPDM ^ ' )■ 

For the decay signal (right-hand envelope), the ex- 
ternal parts of the halo contribute more significantly to 
the signal because the quantity to integrate is shallower 

(pdm — > '''~^)- We note that in the MCXC catalogue, the 
mass ratio of the most massive galaxy cluster to the least 
massive one is ~ 100, so that these clusters have different 
concentrations i?vir/^sj hence a diff'erent rescaled angular 
size avir/c«s- It does not impact on the DM annihilation 
signal from the smooth halo (because this signal is strongly 
centrally concentrated), but it leads to a small spread for a 
decay-like signal. We checked that using an Einasto profile 
(Springel ct al. 2Q0i) instead of a NEW for the DM halos 
does not change the curves significantly. 

For CRs (black envelope), the trend is similar to the 
decay case because Pgag — > j.-i05^ rpj^^ ^-^^^^ between the 
CRs and decay envelopes comes from the gas being more 
centrally concentrated than DM. 
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Fig. 1. Signal fraction fxiv) (for process X) as function 
of the reduced integration angle ?; = a-mt/ois- For the sake 
of presentation, only the envelopes for the 1743 clusters of 
the catalogue are shown. For annihilation, the signal for the 
smooth only is also shown (the substructure contribution 
only is similar to the decay case). 



The DM annihilation from substructures is directly pro- 
portional to the integration of the substructure spatial dis- 
tribution dNsuhs/dV, which we have chosen in this work 
to follow the smooth DM profile. Because of this choice, 
the signal from substructures follows the same dependence 
as DM decay and is therefore not repeated in the figure. 
Similar to the smooth annihilation signal, changing this dis- 
tribution does not impact on the conclusions significantly. 
As underlined, the more complex (i.e. non-universal) scal- 
ing of the total annihilation signal shown in Fig. [T] comes 
from the fact that both the smooth halo and the sub- 
structures contribute to the signal. As shown in Fig. 4 of 
iNezri et al.l(l2012[ ). the distribution of boosts for the MCXC 
sample depends on the substructure configuration, and it 
can vary from 1 (no boost) to about 100 for a Phoenix- 
like configuration. In the former case, the angular depen- 
dence would be close to the smooth only annihilation case, 
whereas for the latter case, it would show a decay-like an- 
gular dependence. 

To conclude, whereas the angular dependence is some- 
times argued to be a good diagnosis to disentangle a CR 
signal from a DM signal, we draw attention to the fact 
that the situation is not that clear cut for galaxy clusters. 
Indeed, a CR signal is more centrally concentrated than 
a DM decay signal. However, for the annihilation case, it 
could be more concentrated than the CR signal, close to 
the DM decay case, or in-between. 



3.2. Number of objects per flux bin 

We now turn to the behaviour of the number of objects 
found given their flux, i.e. the logiV — logi^ dependence, 
where F oc Ax is the actual observed flux (it is referred to 
as the log^ — ]og Ay dependence below). Let us deflne the 
relative flux A^^ of the i — th galaxy cluster, i.e. the flux 
relative to that of the brightest object of the catalogue for 




Fig. 2. Number of galaxy clusters per logarithmic bin of 
Ax (sec Eq. Q), for the CR component (solid black), DM 
annihilation (dotted red), and DM decay (dashed blue). 
The four panels correspond to four integration angles aint 
(from top left to bottom right, 0.01°, 0.05°, 0.1°, and 0.5°). 
The index s of the power law N ex A^ flt is shown for 
ttint = 0.1° for the three components. 



the process X at stake: 



A\ 



A 



X 



X 



brightest 



A 



(7) 



X 



Such a definition allows us to visually compare the be- 
haviour for the different X processes. 

This is shown in Fig. [5J where the number of ob- 
jects N per logarithmic bin (of Ax) is plotted for the 
MCXC_sample. At first glance, the value of the slope s 
{N ex Al^^ oc Ax'^) in the region beyond the distribution 
peak can be used as a diagnosis of the signal origirQ. We 
observe that the behaviour of the DM annihilation and CR 
components is very similar. This is tentatively explained 
by the fact that the CR emissivity Eq. (g]) is oc p^^s ^ ^^^ 
and DM annihilation is oc Pdj^. For the latter, the slope 
Sannihii. varics in the range [sdocay/2, Sdocay]: for small inte- 
gration angles (top left), £'dccay oc poM and £annihii oc p^j^ 
and Sannihii. ~ Sdocay/2 (substructurcs scarccly play a role 
in the central regions of the cluster), whereas for inte- 
gration angles encompassing the whole cluster (bottom 
right), the substructure contribution becomes dominant 
and Sannihii. ~ Sdocay Actually, the annihilation signal at 
any given integration angle is sensitive to the DM substruc- 
ture configuratio n (as explained in the previous section, see 
also Sect. 3.4 of iNezri et al.l[2012f ). so that the values of 
Sannihii. sliowu in Fig. [5] are very uncertain (they are brack- 
eted by the value obtained for the CR case and the value 
for the decay case). 



^ The catalogue is not complete across the range of masses 
it samples, explaining the peaks seen in Fig. [2] (the selection 
function is undetermined, beyond the scope of this paper, and 
not accounted for in the discussion). 
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Fig. 3. Cumulative of the signal above a threshold value Athrcsh for the signal-limited regime (left panel, corresponding to 
Eq. ([H) and for the background- limited regime (right panel, corresponding to Eq. ([9])). The three line -widths correspond 
to three integration angles, and the three styles and colours correspond to CRs, DM annihilation, and decay. An arbitrary 
normalisation is set in order to ease the comparison by eye of the different signals. 



The actual value of s is difficult to explain from simple 
arguments because it comes from the interplay of the pro- 
file, normalisation, and integration angle for each cluster. 
This analysis nevertheless sho-ws that a CR origin can be 
disentangled from a DM origin in the case of a decaying 
DM candidate. 



3.3. Stacked signal 

As deduced from Fig. [31 the ratio bet-ween the brightest 
and next-to-brightest clusters could also be an indicator 
of the 7-ray signal origin (it is the largest for CR origin 
and the lo-west for decay at 0.5°, bottom right-hand panel). 
This is, however, sensitive to the individual modelling of 
the clusters, especially for the annihilation case because of 
the extra uncertainty from DM substructures. In that re- 
spect, this ratio criterion (or also using the ratio of the CR 
signal to the DM annihilation signal, see lPinzke et al.l[201l[ ) 
should be used -with caution. The cumulativ^H of the sig- 
nal also suffers from similar uncertainties. Ho-wever, as the 
number of stacked objects increases, the sensitivity to the 
modelling of a given halo becomes less crucial. The possibil- 
ity to identify the origin X of a signal from the behaviour 
of the stacked signal is very similar in spirit to comparing 
the slope s of the logiV — log^x relation. Ho-wevcr, if the 
signal is only seen by means of stacking, it may not be pos- 
sible to get s, but still possible to say something about the 
origin from this detection. 

The promise of a stacking analysis for a signal X de- 
pends on the slope s of the logA^ — log^x relation. If 
the number of objects increases faster than the signal de- 
creases, a better sensitivity should be reached compared 
to the single-source analysis case. This w as shown to be 
the case fo r decay in Combet et al.l (|2012h , but not for an- 
nihilation (iNezri et al.ll2012l ). A detailed discussion of de- 
tectability in the light of the Fermi-LAT and CTA-like in- 
struments was also presented in that paper. Here, we extend 
the analysis to the CR case, and comment on the differences 



between the various signals. The simple approach used in 
ICombet et al.l ()2012[ ) is enough to capture the trends and 
draw some conclusions, i.e. we consider the generic extreme 
case of a signal- and a background-limited observatior(f|. 

As in the previous section, we work with the normalised 
flux Ax defined by Eq. ([7|) to ease the comparison between 
the different X origins. The cumulative for a signal-limited 
observation (left-hand panel) corresponds to 

^Sig-lim(aint) OC 2J ^x(aint), (8) 

while that for a background-limited observation is 



'CBkgd-lim(aint) « 



/Csi 



Sig— liin 



(ctint) 



aintV^(>A„,_0 



(9) 



where -/V(>^ - ) ^^ ^^^ number of clusters satisfying A^-^{= 
^^/^^"s''*°'*) > ithrcsh- This is shown in Fig. |3] for the 
MCXC cataloguco- For display purposes, the curves for var- 
ious processes X are normalised to the asymptotic value 
(when all objects are stacked), whereas the curves for dif- 
ferent integration angles aint are shifted from one another. 
Regardless of the integration angle and whether the 
regime is signal or background limited, stacking leads to 
a significantly larger increase of the signal for DM decay 
than for the CR case. Indeed, as shown in ICombet et al.l 
(|2012l ) for the DM decay component, improvement with re- 
spect to the brightest source can reach at best (with the 
MCXC catalogue) a factor of ^ 100 for background-free 
instruments, but a factor -^ 5 is expected for more realistic 



'' The cumulative is defined to be the cumulated flux of all 
objects brighter than a threshold flux 



* The detectability of galaxy clusters depends both on the in- 
tegrated signal and on the much higher level of charged particle 
and diffuse astrophysical 7-ray backgrounds. In the background- 
limited regime the best approach is to maximise the signal-to- 
noise ratio. For a uniform background, it corresponds to the 
cumulative of Ax divided by the square root of the background 
which is OC af^.^^N. 

® The plateau reached by the cumulative signal on the left- 
hand panel is due to the drop of the number of objects seen in 
Fig. El 
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background-limited instruments. For both the CR and DM 
annihilation component, the improvement is respectively 
^ 50 (background free) and ^ 2 (background limited) for 
c^int ^ 0.1°. A detailed analysis based on a realistic in- 
strumental response for the Fermi-LAT and CTA instru- 
ments shows that this factor of two (for annihilation) is 
only achieved for Fermi -LAT, there is no improvement with 
CTA-like instruments (jNezri et al.l 12013) . Hence a similar 
result is expected for the CR-induced origin. Although in 
principle the behaviour of stacked signals at different an- 
gles for the DM annihilation case can be disentangled from 
the two other cases for the signal-limited regime, the uncer- 
tainties on the substructure distribution are such that the 
spread generated on the stacked annihilation signal is en- 
compassed by the CR and DM decay case and thus cannot 
be distinguished for now. 

4. Conclusion 

The first two papers of this series investigated the potential 
benefit of a stacking analysis of clusters of galaxies for DM 
detection. Relying on the large st uniformed X-ray ca talogue 
to date, the MCXC catalogue (jPiffaretti et al.ll201lh . it was 
shown t hat a DM decay sig nal benefits from a stacking 
analysis (|Combet et al.ll2012|) , but th at the situation is less 
promising for the annihilation signal (jNezri et al.ll2013 ). In 
this third and last paper of the series, we extended the 
analysis to the expected CR-induced component, hoping to 
put to the forth a difference that may be used as a novel 
diagnosis to identify the origin of the signal when detection 
in such objects becomes routinely available. 

Based on the slope s of the log N — log Ax behaviour 
(or equivalently on the benefit gained by stacking the sig- 
nal from many objects), we have found that a DM origin 
can be identified against an astrophysical background in 
the case of decaying DM, but not for annihilating DM. 
Though, in principle, combining the angular-dependence 
information and the stacked signal at different integration 
angles gives clues as to the signal origin, the existing un- 
certainties on the DM substructures (which may contribute 
significantly to the annihilation signal) prevents us from 
reaching any clear conclusion. Furthermore, whereas it has 
been sometimes argued that the angular dependence can 
be used to disentangle DM decay from DM annihilation, we 
have shown for galaxy clusters that i) there is not strictly 
a universal angular dependence for the annihilating signal, 
and ii) that the very uncertainties on the DM substructures 
prevent us from predicting exactly what this dependence is. 
Hence, decaying DM is the simplest scenario to test and/or 
exclude, because it is the most sensitive to stacking and the 
angular dependence of several objects can also be stacked 
(with a rescaling). 

The robustness of this analysis would benefit from a 
more thorough investigation that takes into account vari- 
ous sources of uncertainties, such as the DM and gas pro- 
file modelling. For instance, the intrinsic scatter observed 
in gas density profiles at the centre of clusters (due to gas 
physics, interactions with the central galaxy and with the 
often present active nucleus, the dynamical state of the clus- 
ter, etc.) could affect the CR signal signature. Similarly a 
more systematic check of the dependence of the DM anni- 
hilation signal on the clump parameters (for the substruc- 
tures' contribution) is certainly needed. The calculations 
presented here have been done in the idealised scenario 



where all clusters have a dominant signal component of the 
same nature, whether it be DM decay, annihilation or astro- 
physical. The effect may be even less clear if this is not the 
case. Nevertheless, this may be a useful test that will com- 
plement those already suggested in the literature (spectral 
feature, spatial-dependence of the signal, multi- wavelength 
analysis, angular power spectrum). 
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